Adiabatic creation of coherent superposition states via multiple intermediate states 
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We consider an adiabatic population transfer process that resembles the well established stimu- 
lated Raman adiabatic passage (STIRAP). In our system, the states have nonzero angular momen- 
tums J, therefore, the coupling laser fields induce transitions among the magnetic sublevels of the 
states. In particular, we discuss the possibility of creating coherent superposition states in a system 
with coupling pattern J = <4> J = 1 and J = 1 J — 2. Initially, the system is in the J = state. 
We show that by two delayed, overlapping laser pulses it is possible to create any final superposition 
state of the magnetic sublevels |2, —2), |2, 0), |2, +2). Moreover, we find that the relative phases of 
the applied pulses influence not only the phases of the final superposition state but the probability 
amplitudes as well. We show that if we fix the shape and the time-delay between the pulses, the 
final state space can be entirely covered by varying the polarizations and relative phases of the two 
pulses. Performing numerical simulations we find that our transfer process is nearly adiabatic for 
the whole parameter set. 

PACS numbers: 42.50.Hz, 42.65.Dr, 32.80.Bx 
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I. INTRODUCTION 

Adiabatic methods in quantum mechanics play a fun- 
damental role in the treatment of level crossing problems. 
The famous Landau-Zener system has been followed 
by several other models which were defined in two-level 
systems 0, S 0, 0! ■ These models induced many applica- 
tions in atomic and molecular physics, for recent review 
seeRef. @. 

It is common in all of the adiabatic systems that 
the Hamiltonian is time-dependent. The instantaneous 
eigenstates define an adiabatic basis. Suppose that at the 
initial time, one of the adiabatic states coincides with the 
initial state of the system. In most cases the adiabatic 
state that belongs to the eigenvalue zero is utilized for the 
transfer. If the system varies with time slowly enough, 
then its state vector follows adiabatically the adiabatic 
eigenstate and it will go over smoothly to the desired fi- 
nal state. The adiabaticity of the process guarantees its 
robustness with respect to fluctuations in the parameters 
of the couplings. In the limit of ideal adiabatic transfer, 
the other eigenstates are not invoked in the evolution 
and their population remains zero throughout the whole 
time. The non-adiabatic corrections tend to involve these 
states, but the magnitudes of such couplings are supposed 
to be small compared with their energy separation from 
the adiabatic states. 

During the last decade, a promising new population 
transfer technique was introduced: the stimulated Ra- 
man adiabatic passage (STIRAP) process realizes effi- 
cient population transfer in three- level systems, see Ref. 
and references therein. Moreover, it turned out that 
this scheme can be applied successfully to the problem of 
manipulating and creating coherent state superpositions. 



Such superpositions are the desired initial states for many 
modern quantum applications including information pro- 
cessing and communication. The original STIRAP pro- 
cess has thus been utilized to create coherent superposi- 
tions in three- and four- level systems |al2lllfllllllll3l an d 
to prepare A-component maximally coherent superposi- 
tion states 13]. The four-level (tripod) STIRAP estab- 
lishes such applications as qubit rotation |14| and density- 
matrix measurement |15). 

Recently, a multilevel STIRAP scheme has been pro- 
posed [llJ 17] to create A— level coherent superposi- 
tion states. This model is particularly interesting, since 
the Hamiltonian has several adiabatic eigenstates which 
belong to the eigenvalue zero. The diabatic couplings 
among these eigenstates influence significantly the dy- 
namics of the system. However, the eigenstates in the 
zero eigenvalue subspacc arc decoupled from the eigen- 
states belonging to nonzero eigenvalues provided the 
time-evolution is adiabatic. In Ref. |lfj( the Optimal 
Control Theory has been applied to find the optimal 
pulse sequences which create a prescribed final super- 
position state. It has also been verified that the transfer 
process occurs in the zero-eigenvalue subspace. In a sub- 
sequent work |17| an analytic solution has been presented 
for a five-level system. 

When the STIRAP concept is realized in atoms or 
molecules one encounters that the quantum states with 
angular momentum J have magnetic sublevels and these 
sublevels are coupled by the applied laser pulses, and a 
rather complicated coupling-pattern may occur. This sit- 
uation has been studied in Refs. [l8L ll^. l2Cj. In these 
papers the authors have discussed the possibility of trans- 
ferring population from a well defined initial state to a 
single but arbitrarily chosen magnetic sublevel. 

An other realistic problem is the presence of multiple 
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or continuum intermediate states [2III2II2II2 I I2 I I2I 
I2H |28| . In Ref. [2l| the continuum has been modeled by 
infinite number of equidistant discrete levels which were 
coupled by unique couplings to a single initial and final 
state. Later, it has been shown that the previous model is 
too simple, the physics is more complicated in a realistic 
situation |23 • A detailed study of the discrete model has 
been presented in Ref. [23J. It has also been shown that 
efficient po pula tion transfer is indeed possible through a 
continuum [H |H |H El E| . 

In this paper we will study a multilevel A scheme simi- 
lar to the one discussed in Refs. [Hi. IT^ . I20I ] . However, in 
contrast with the previous works where the population 
transfer took place between a single initial and a single 
final state, our aim is to create a multilevel final super- 
position state. In our model, the initial state has angular 
momentum J=0, the intermediate state has J=l and the 
final superposition state is formed from the magnetic sub- 
levels of a state J=2. We will design a STIRAP-like pro- 
cess, i.e. the population transfer is realized by delayed, 
overlapping laser pulses. We define the Hamiltonian cor- 
responding to the process, separate the adiabatic states 
belonging to the eigenvalue zero and find the conditions 
for the successful population transfer. We will study how 
does the final state depend on the relative phases of the 
applied laser pulses. The robustness of the transfer pro- 
cess will also be discussed. 

The organization of the paper is as follows: In Section 
ITTl the physical model of our six-level STIRAP is pre- 
sented and the corresponding Hamiltonian is introduced. 
The Hilbert space is separated to a zero-eigenvalue sub- 
space and to a nonzero-eigenvalue subspace. The eigen- 
states belonging to the zero eigenvalue are determined 
explicitly In Section Hill we study the possibility of pop- 
ulation transfer by means of delayed laser pulses. We 
show numerically that nearly adiabatic evolution is pos- 
sible. We also discuss the robustness of the transfer pro- 
cess. We summarize the results in Section llVl 



II. SIX-LEVEL SYSTEM 

Let us consider the six-level A configuration shown in 
Fig. ^ The system consists of the magnetic sublevels 
of states with total angular momentum J=0, J—l and 
,7=2. The ,7=1 states have higher energy than the J—Q 
and J=2 states. Initially, only the J=0 state is popu- 
lated. Our goal is to transfer this population to the sub- 
levels of the J=2 state and to obtain a prescribed final 
superposition state. We are going to design a STIRAP- 
like population transfer process. To this end, we apply 
two sets of laser pulses, the corresponding Rabi frequen- 
cies are denoted by Sls± and Q p ±- The subscripts S and 
p refer to "Stokes" and "pulse", while the indices +/— 
correspond to right- and left-hand circular polarizations, 
respectively. Obviously, in this configuration only the 
sublevels with m— — 2, 0, +2 can be populated in the fi- 
nal state. This linkage can be realized experimentally in 
Neon atoms, see Ref. [2(j. The more general case, where 



all magnetic sublevels of the J=2 state are included into 
the final superposition state will be discussed elsewhere. 
The frequency and the phase of the laser pulses are kept 
constant during the interaction, only their amplitudes are 
time-dependent . 

In our model we have to determine the Rabi frequencies 
in an angular-momentum basis |29|. According to the 
Wigner-Eckart theorem, the general formula for the Rabi 
frequency between two levels of angular momenta J\ and 
J 2 and magnetic quantum numbers M\ and M 2 is given 
by 
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where the electric field strength is £ , the reduced dipole 
matrix element is given by ( Ji \ \d\ \J 2 ), e.s axe the irre- 
ducible tensor operator components of the polarization 
vector of the field, and [. . . ] denotes the 3 — j symbol. 
We apply left- and right-hand circularly polarized light, 
therefore, in the summation s = ±1 and e±\ = 1. 

The two pump pulses couple the level 1 00) to the levels 
|1,±1). The corresponding Rabi frequencies can be ex- 
pressed as £l p ± = consti x £*± . The right-hand circularly 
polarized Stokes pulse couple the levels |1, — 1) |2, — 2) 
and |1, +1) |2, 0). For the Rabi frequencies we obtain 

f2s + (l, — 1; 2, — 2) = q x const 2 x £g + , (2a) 
O s+ (l,+l;2,0) = const 2 x£| + , (2b) 

where q is the ratio of the corresponding 3 — j symbols. 
In our case q — y/d. Similarly, the left-hand circularly 
polarized Stokes pulse couple the levels 1) 1 2 , 0) 
and |1, +1) |2,+2), and the Rabi frequencies are given 
by 



fi lS _(l J -l;2,0) 
fi s _(l,+l;2,+2) 



const2 x £g_ , (3a) 
q x const2 x £g_ . (3b) 



In summary, the Rabi frequencies can be expressed as 



Q p+ (0,0;1,-1) 


= A expitpA , 


(4a) 


fi„_(0,0;l,+l) 


= B exp iipB , 


(4b) 


n s+ (l,-l;2,-2) 


= qC exp iepe , 


(4c) 


n s+ (l,+l;2,0) 


= C exp iipe , 


(4d) 


n s _(l,-l;2,0) 


= D exp iifD , 


(4e) 


n s _(l,+l;2,+2) 


= q D expiipz, . 


(4f) 



The phases ipx are constants, and the envelope func- 
tions X are non-negative. The laser pulses may be de- 
tuned from resonance by a frequency A, but we require 
that they are at multi-photon resonance. The interaction 
Hamiltonian in the basis { |0, 0) , |1,-1), |1,+1), |2, —2), 
|2,0), |2,+2)} reads 
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H 



A exp tip a 

Bexp—iipB 

qC exp —iipc 

D exp —iifD 











B expi(fiB 
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C exp —iipc 
qD exp —itfiD 





qCexpitpc Dexpipjj 

C exp z</?c ?-D exp i</>d 





(5) 















The time-evolution of the system is governed by the 
Schrodinger equation 



where 



(6) 



In a STIRAP process, the Hilbert space is split to 
dark and bright subspaces and the system evolves in the 
dark subspace. These subspaces are formed of the in- 
stantaneous eigenstates of the time-dependent Hamilto- 
nian. The dark states belong to the eigenvalue zero and 
they do not have projection to the excited state(s) [30| . 
The bright states belong to nonzero eigenvalues and they 
overlap with the excited state(s). In the adiabatic limit, 
the system evolves entirely in the dark subspace. Now 
we continue the analysis of the Hamiltonian Eq. (J5J in 
this direction. First, we determine the dark and bright 
subspaces of the Hilbert space. By applying a diagonal 
unitary transformation Ui to the system, the number of 
independent phases in the Hamiltonian can be reduced 
to one: 
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where 



U x = diag[e- j(vA+vc) ,e-' iVC ,e^ D - 2 ^ c) ,l, 



e i(<pn-<pc) e i(2yD-2(pc)l 



and the remaining phase factor 

if = ip B - p A + ip c - <p D 



(7) 



(8) 



(9) 



It will turn out that this phase factor is indeed relevant, 
since it affects the probability amplitudes of the final 
state. 

As a next step we are going to find a time-dependent 
unitary transformation that defines a new basis, in which 
the dark and bright subspaces of the Hilbert space are 
separated. By applying rotations and diagonal unitary 
transformations in a well-chosen order we arrive at the 
following unitary transformation: 



u = U 4 U 3 U 2 , 



(10) 
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(11) 



and 



U3 = diag[expi£, 1, 1, 1, expi£, 1] , 



(12) 



U 4 = 



cos a 








1 
1 
1 



sin a 






- sin (3 sin a cos [3 sin (j cos a 
_ — cos (3 sin a — sin f3 cos (3 cos a 



where we made use of the definitions 



k 
x 

y 

tan a 
tan f3 

c 



C/D, 



= V^ + fcV + i/V^ + i, 

= -q/Vk 4 + l, 

(fc 4 + l)yD 

^A 2 k 2 + B 2 - 2ABk cos ip ' 

_y j A 2 +B 2 k e + 2ABk 3 cosy ^ A 

- -- C0 W A 2 k 2 + B 2 -2ABk cos V { ' 



(13) 

(14a) 
(14b) 
(14c) 

(14d) 



arctan 



B sin ip 



£ + arctan 



Ak — B cos <p, 
Bk 3 sin ip 



A + Bk 3 cos ip 



(14f) 
(14g) 



If we transform the Hamiltonian Eq. Q with the uni- 
tary transformation U in Eq. (|1C1|) we obtain 



Hx = UHxU 1 



(15) 



which has two rows and two columns with zero elements. 
This means that the unitary transformation U defines 
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two dark states. The first one is given by the fourth row 
of U 2 





D 2 

-qCD 
C 2 



\Di) = 



VC 4 + L> 4 + q 2 C 2 D 2 



(16) 



The one-dimensional subspace corresponding to that 
dark state is not altered by the other unitary transfor- 
mations C/3 and U4, and it is a superposition of the levels 
with J=2. The sixth row of the matrix U defines the 
second dark state 

— cos (3 sin a exp j£ 





\Do) = 



■— sin j3 exp i£ 



m sin 3 exp iC + cos /3 cos a 

qx q 



2 .3 

- ^^r- sin /3 exp iQ — | cos /? cos a 



(17) 



The states and I-D2) are indeed dark states, since 
they do not have components among the J—l levels. 

To complete the study of the dark and bright sub- 
spaces of our system we should diagonalize completely 
the Hamiltonian Hi in Eq. 115fl . This leads to rather in- 
volved formulae, however, for our purposes it is enough 
to show that the eigenvalues of the 4x4 nonzero block 
of Hi are nonzero and the corresponding eigenstates are 
bright states. We have verified numerically using Maple 
[3l| that the previous assumptions are true. We conclude 
that the states given by Eqs. I|16fl and (|17fl are the two 
dark states of our system and it has four bright states. 
The dark subspace is degenerate, so particular attention 
should be payed when one studies the time-evolution of 
the system [3, 0, 0- In general, for time-dependent 
dark states there appear a nonadiabatic coupling term 
in the transformed Hamiltonian which couple the dark 
states. We will return to this issue in the Sec. IIIII We 
can transform the Schrodinger equation 10 to the time- 
dependent basis defined by Ui and U in Eqs. © and ltTU|) 



through the relations 



H = Hi + iUU^ , 



(18) 



(19) 



where Hi is given by Eq. I|15fl ■ In this form the dark and 
bright subspaces are separated. Since we want to design 
a STIRAP-likc population transfer process, in the follow- 
ing we focus our attention to the dark states. The bright 
states are needed when one wants to study the coupling 
between the dark and bright subspaces, but since the for- 
mulae for the bright states are too complicated, we shall 
do it in an other way. 



III. ADIABATIC POPULATION TRANSFER 

Now we make a further assumption: We assume that 
the ratios A/B and C/D are constant during the time- 
evolution, i.e. the envelope function of the two pump 
pulses are the same, only their amplitudes and phases are 
different. Similarly for the Stokes pulses. Consequently, 
the parameter k defined in Eq. (|14afl is constant. There- 
fore, the dark state |Di) is constant throughout the whole 
time. This simplifies considerably the further analysis. It 
follows immediately, that though the dark subspace is de- 
generate, the two dark states do not mix with each-other 
during the interaction, since (-D2I-D1) = 0. Therefore, the 
Hamiltonian H has a 2 x 2 zero matrix block acting on 
the dark subspace. 

The dark state I-D2) has components among the J=0 
and J— 2 states as well. Hence, it is a good candidate 
to realize the transfer process under consideration. Ini- 
tially the system is in the state |0,0). For implementing 
a STIRAP process, that state should be a dark state of 
the system. Since the state \Di) cannot participate in a 
superposition forming the initial state, we require ID2) 
to coincide with the state |0, 0) at early times. This con- 
dition is satisfied if C 2 + D 2 ^> A 2 + B 2 , which means 
that the pulses C and D arrive before A and B. Though 
we have not presented explicitly the bright states of our 
system, it is sure that they have no projection to the ini- 
tial state, since they are always orthogonal to the dark 
states and \D 2 ) = |0, 0) initially, if the above condition is 
fulfilled. This argument ensures that the initial state of 
the system lays in the dark subspace which is necessary 
for a STIRAP process. 

We have required that the final state of the system 
is a linear combination of the levels with J=2 and 
m= — 2, 0, +2. The dark state I-D2) has these compo- 
nents as well, so the previous condition is true only if 
the first element of the state I-D2) is zero at the end of 
the adiabatic evolution. Consequently, at late times it is 
eligible that C 2 + D 2 < A 2 + B 2 , therefore, the pulses 
C and D terminate before A and B. 

The previous conditions suggest that for successful 
population transfer we need counterintuitive time order- 
ing of the pump and Stokes pulses, similarly to the orig- 
inal STIRAP. Now we consider the possibility of adia- 
batic population transfer using the connectivity argu- 
ment of Ref. 0| : The system Hamiltonian has a con- 
stant zero eigenvalue which is degenerate. However, as 
we have shown in the beginning of this Section, the nona- 
diabatic coupling between the two degenerate dark states 
vanishes because one of them is constant. As we have 
seen in the previous paragraphs, the time dependent dark 
state coincides with the initial state of the system at the 
beginning of the process when only the Stokes pulse is 
present. Similarly, this dark state lays in the required 
final state subspace at the end of the process when only 
the pump pulse is present. Therefore, the connectivity 
requirements are satisfied |19| . we have got a dark state 
that connects smoothly the initial state of the system 
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with the final state. However, besides these conditions 
we have to fulfill adiabaticity as well. We will discuss 
this issue below. 

The considered six-level system can be regarded as four 
coupled three-level systems: For A = and C = the 
pulses B and D define a three-level ST1RAP process be- 
tween the levels |0,0), |1,+1) and |2,0). For A = and 
D = the pulses B and C define a three-level ST1RAP 
process between the levels |0,0), |1, +1) and |2, +2). For 
B = 0, C = and B = 0, D = the situation is similar. 
The envelopes A, B, C, D can be parametrized in a form 
that emphasizes the presented qualitative picture: 



.4 
B 

C 
D 



R p cos 77 , 
R p sin x] , 
Rs 

; I 

Rs 



(20a) 
(20b) 

(20c) 
(20d) 



where the constant angles 77 and v are in the interval 
[0,7r/2]. These angles define the polarization of the laser 
pulses. The six-level system simplifies to a three-level A 
system when 77, v = 0, w/2. For arbitrary constant an- 
gles 77 and v fixed couplings exist between the individual 
three-level STIRAP processes. 

If the time-evolution is adiabatic, then the dark and 
bright subspaces do not mix during the whole interaction 
time |9tll6ul7|. Therefore, if the initial state of the sys- 
tem lays in the dark subspace, it will remain there. The 
physical consequence of the adiabatic evolution is that 
the excited states |1,±1) are only minimally populated 
throughout the transfer process. As we saw in the previ- 
ous paragraph, when 77, f = 0, tt/2 the system reduces to 
the ordinary three-level STIRAP. In this case the adia- 
baticity constraints are simple: Let us take 77, v = 0, 7r/2, 
then adiabaticity requires 



. 2 
sin g 



cos g 



cos 2 g 



sm g 



<ifi, (21) 



where f2 



ng, tan?9 = Clp/ntts, ( K = !><?)) 
and tan2g = f2/A. These conditions say that we need 
large pulse amplitudes and smooth, long pulses. For 
a complete study we would need the explicit form of 
the bright states \Bk), k = 1 . . .4 of the Hamiltonian in 
Eq. 10. Then, adiabaticity requires to satisfy the follow- 
ing four equations: 



test adiabaticity numerically for a certain choice of the 
pulses. In the numerical simulations presented below we 
used these equations to test adiabaticity. Besides these 
conditions, below we show that a simplified approach is 
also possible: we formulate a sufficient condition for adi- 
abatic evolution in the six-level system. 

Our physical intuition suggests that if the three-level 
STIRAP processes for the special polarizations 77, v = 
0, 7r/2 are adiabatic, then the six-level STIRAP process 
will also be adiabatic for constant angles 77 and v. This is 
a sufficient condition for adiabaticity. It is independent 
of the polarization of the pulses, in fact, we apply the 
adiabaticity condition of the original three- level STIRAP 
to our six-level scheme. We have performed numerical 
simulations to reveal whether this assumption is valid or 
not. 

In Figs. |2| and the numerically calculated time evolu- 
tion of the populations are shown for two different param- 
eter sets. The pulses are Gaussian. The two parameter 
sets differ only in the phase tp. The final populations, de- 
noted by /1, fi and f% in the figures, are different, hence 
the phase ip has indeed impact on the magnitude of the 
probability amplitudes in the final state. We interpret 
this phase dependence as a manifestation of quantum in- 
terference among the rivaling coupled population transfer 
processes. 

To test adiabaticity in the previous examples, we eval- 
uated Eqs. (|22|l for the parameters of Fig- El The results 
are shown in Fig. 0] It can be seen that for that time- 
interval when adiabaticity is expected the conditions are 
well fulfilled. The non-adiabaticity of our six-level STI- 
RAP process can be also quantified by the maximum 
population of the bright subspace during time-evolution. 
This non-adiabaticity is 0.014 and 0.011 for the two pre- 
sented processes. The small values show that the mix- 
ing between the dark and bright subspaces is negligible, 
therefore, the time-evolution of the system is adiabatic. 
The efficiency of the population transfer can be defined 
as the sum of the norm of the final states with J = 2. For 
the two presented processes it is found to be 99.98% and 
99.97%, showing that it is possible to perform a STIRAP 
process in the considered six-level system. 

The final state of the system in the adiabatic limit 
can be obtained analytically from the expressions of the 
unitary transformations U and U± at the beginning and 
at the end of the population transfer 



\(D 2 \B k )\ < |e*-eo| 



k = 1...4, 



(22) 



where s k is the adiabatic eigenenergy associated with 
\Bk), and £0 = belongs to the dark state \D 2 ). These 
equations yield the necessary and sufficient conditions 
for the pulses to be fulfilled in order that the evolution 
satisfy adiabaticity. Although the bright states \B k ) are 
not available analytically, these equations can be used to 



\^ f ) = ulu]u i u 1 \^ i ) 



(23) 



provided that lays in the dark subspace. In our case 
the initial state of the system is \ipi) — |0, 0). After some 
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calculation we obtain 



Q 
Q 


(_J0si n/ 3 e -< + Ul^cos^e^^c-a 

K - y C0S l3) e -i('PA-'Pa+2 V n-S) 



sin /3e 

V qx 1 



(24) 

where the parameters k, x, y, f3, £, £ are defined by Eqs. 
(I14JI . and the parameters /3, £, and £ are evaluated at the 
final time. The equations Il2i'>[l and H24[) together with 
(|14|l show that our population transfer process is robust 
with respect to small fluctuations of the experimental pa- 
rameters such as pulse shapes and pulse areas. However, 
the relative phases of the pulses influence not only the 
phases in the final superposition state but the probabil- 
ity amplitudes as well. 

It is necessary to verify that all possible linear combi- 
nations of the states with J~2 and m= — 2, 0, +2 can be 
created by the considered STIRAP process. Suppose, we 
have a prescribed final state 



\i>c 



o 




cie l< - 
cae* 
c 3 e 



vp 3 



(25) 



and we have found a pulse set with ipA = 0, <ps — <p, 
ipc = 0, ipo = for which the final populations of the 
state \ipf) are the same as in the prescribed state, and 
the equation 



<y5i + ¥3 - 2<y3 2 = arg|^/)4 + arg|-0/)e - 
— 2arg|^/)5 mod 2tt 

holds. By adjusting the phases of the pulses to 

ip A = arg|-0/)5 - <^2, 

ip B = <ys + arg|V>/)e - ¥>3, 

ipc = arg|-0/)5 - arg|V'/)e + <^3 - ^2, 

fD = 0, 



(26) 



(27a) 
(27b) 
(27c) 
(27d) 



the phases of the final state change to the phases of the 
prescribed state. The phase ip defined by Eq. |0 is invari- 
ant under the presented phase adjustment, so the final 
populations remain the same. It follows that for deciding 
whether all prescribed final states can be reached, it is 
enough to show that all possible amplitudes (ci, 02,03) 
and all possible values of 



+ ¥>3 — mod 2tt, 



(28) 



can be reached by pulse sets with tps = <p and ipA — 
<Pc = ¥>d = 0. 

The amplitudes (c\, C2, C3) satisfy of + c\ + c\ = 1. This 
triplet defines a point on the one-eights of the surface of 



the unit sphere, since the restriction ci, 02,03 > must 
be satisfied. The surface points can be parametrized by 
two polar angles 9 £ [0, tt/2] and \ G [0,tt/2], therefore 
we have 



c\ = cos 9 , 

02 = sin 9 cos x ■ 

03 = sin 9 sin \ ■ 



(29a) 
(29b) 
(29c) 

If we fix the shape of the pulses and the time-delay be- 
tween them, then we have three independent parameters 
which can be adjusted freely: the polarizations 77 and 
v and the combination ip of their relative phases. The 
final state can be characterized by the triplet [9, x, S). 
Mathematically, our population transfer process assigns 
to each triplet (77, v, ip) a triplet {9, \, S). If these points 
fill the cube [0 ... tt/2, ... tt/2, ... 2w], then all pre- 
scribed states can be reached by choosing a correspond- 
ing triplet (rj, v, ip) and performing the transformation 
defined in Eq. I|27|) on the initial phases of the pulses. 

We have performed numerical simulations in order to 
answer the above question. The angles 77, v and ip 
were incremented by 0.022 in their domain range. We 
have used the envelope functions R p — 15exp[— ((f — 
1.8)/2.82) 2 ] and R s = 150 + q 2 exp[-((* + 1.8)/2.82) 2 ] 
for the pump and Stokes pulses, respectively. We have 
verified that these pulses fulfill the adiabaticity condi- 
tions Eq. I|21|) for the field polarizations r\,v = 0, tt/2. 
The resulting triplet set was divided into slices of width 
0.1 in S. Each slice was checked whether it is filled with 
points. Some of the slices are shown in Figs. El CI El 
El El We have found that each slice is equally and quite 
uniformly filled with points. 

We have compared the numerical results with the ana- 
lytical solution given by Eq. (|23|l . Note that the analyt- 
ical solution is valid in the adiabatic limit. In this case 
the increments of the parameters 77, v and ip were 0.015, 
giving a more detailed picture about the distribution of 
the resulting points (9, x, S). The differences between the 
numerical simulation and analytical solution were small, 
below one percent, showing that the chosen pulse-sets re- 
alized a nearly adiabatic evolution. In Fig. 0wc show a 
slice of the distribution for the same parameter set as in 
Fig.0 

The maximum non-adiabaticity, that we measure as 
the maximum population of the bright subspace, encoun- 
tered during the simulation was 0.006. We have also ver- 
ified that the adiabaticity conditions Eq. 1)221) are fulfilled 
for the whole parameter range. Therefore, our assump- 
tion is indeed true: if the population transfer process 
is adiabatic for some special polarizations of the pulses 
77, v = 0, tt/2, i.e. we have an effective three-level system, 
then it remains adiabatic for arbitrary polarizations, i.e. 
all the six states are coupled. This is a sufficient con- 
dition for the pulses to satisfy adiabaticity. Based on 
the results of the numerical simulation and the analyti- 
cal solution, we conclude that all prescribed states can be 
reached by the presented STIRAP-like population trans- 
fer process and the evolution is adiabatic. 
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IV. SUMMARY 

In this paper we have discussed a STIRAP-like popula- 
tion transfer process in a six-level A system. The initial 
state has angular momentum J = 0, the excited state 
has J = 1 and the final state has J = 2. Our aim is 
to create coherent superposition states on the magnetic 
sublevels of the final state. Similarly to the original STI- 
RAP, the couplings are realized by a pump and a Stokes 
laser pulse, which are elliptically polarized in our case. 
The shape of the pulses and the time-delay between them 
are fixed. However, their polarizations and the relative 
phases are freely adjustable but also fixed throughout 
the whole time. We have defined the Hamiltonian of this 
system in the RWA approximation. We have found that 
by time-independent unitary transformations not all the 
phases can be eliminated from the Hamiltonian, there 
remains one. This phase factor proved to be indeed rele- 
vant, since it influences the probability amplitudes in the 
final superposition state. We have determined the dark 
and bright subspaces of our system and we have given 
explicitly the two dark states. We have shown that the 
desired population transfer can be realized by a coun- 
terintuitive pulse sequence similarly to the original STI- 
RAP. We have studied the conditions for adiabatic evo- 
lution in our system: In the case of special polarizations 
of the pulses the coupling pattern reduces to a simple 
three-level system. The adiabaticity conditions are well 
known for this case. We have shown numerically that 
the system evolves adiabatically for general polarizations 
as well if adiabaticity prevails for the special choice of 
the polarizations. We have compared the results of the 



simulation with the analytically calculated final states in 
the adiabatic limit, and an excellent agreement has been 
found. We have verified numerically that the whole final 
state space can be covered by varying the relative phases 
and the polarizations of the two exciting pulses, while the 
pulse shapes and the time-delay between them arc fixed. 

As we mentioned before, the probability amplitudes in 
the final state depend on the relative phases of the laser 
pulses. This is unusual in case of the STIRAP process, 
up to our knowledge that is the first time when the rel- 
ative phases do matter. The phase dependence results 
from quantum interference, because the coupling pattern 
can be considered as a coupled set of three-level STIRAP 
systems. Therefore, the final state emerges from the ri- 
valing STIRAP processes. 

Finally, we comment briefly the experimental realiza- 
tion of our scheme: It seems easy to find an atomic sys- 
tem where the coupling pattern can be attained. Two 
elliptically polarized pulses are required, however, the 
phase shift between the left- and right-rotating compo- 
nents should be adjustable. Practically this seems the 
only additional experimental effort to be done compared 
with the realization of the traditional three- level STIRAP 
process. 
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FIG. 1: The six-level A system. The levels |2, ±1) and |1,0) 
are not coupled to other levels. The remaining levels form the 
six-level A configuration. Levels with J=0 and J=2 are cou- 
pled through the levels with J=l by the pump (p) and Stokes 
(S) laser pulses with right- and left-hand circular polarization 
(+/—). The corresponding Rabi frequencies are Qs± and 
Q p ±, respectively. The pulses are at two- photon resonance, 
but may be detuned from the J—l levels by a certain detuning 
A. Only the level |0, 0} is populated initially. 
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FIG. 2: The numerically computed time evolution of the pop- 
ulations are shown for a set of Gaussian pulses. The pulse 
shapes are the following: R p = 10exp[— ((t — 1.6)/2.8) 2 ], 
Rs = 10^1 + q 2 exp[-((t+ 1.6)/2.8) 2 ]; the polarizations are 
characterized by the angles n — 0.5, v = 1.04; finally ip = 3.34. 
The excitation is resonant A = 0. Time is measured in arbi- 
trary units. 
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FIG. 3: The numerically computed time evolution of the pop- 
ulations are shown for another set of Gaussian pulses. The 
pulse shapes are the same as in Fig. [5] only the phase <p is 
changed to 4.74. 
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FIG. 4: The adiabaticity of the process shown in Fig. [5] The 
upper plot shows the pulses, in the lower figure the ratios 
\sk - £o]/|(.D2|-Bfc)| , k = 1 . . . 4 are plotted. Note that two 
pairs of the eigenstates and matrix elements coincide, hence 
only two curves can be distinguished. Time and frequency 
are measured in arbitrary units. 
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FIG. 5: The slice with 5 = 0.30±0.05 of the cube filled by the 
points (9, x, S) which result from the numerical solution of the 
Scrodinger equation JSJ . We have fixed the shape of the pulses 
and the time delay between them, only their polarizations and 
relative phases were varied. 





7t/4 jt/2 



FIG. 7: Same as Fig.^Jbut 5 = 3.2. 
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FIG. 8: Same as Fig.Qbut the points are computed from the 
analytical solution Eq. 1241 . 





